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ABSTRACT 

We study the radial and vertical stability of dust grains launched with all 
charge-to-mass ratios at arbitrary distances from rotating planets with com- 
plex magnet ic fields. We show that the aligne d dipole magnetic field model 



analyzed by 



Jontof-Hutter and Hamilton! ( 12012| ) is an excellent approximation 



in most cases, but that fundamentally new physics arises with the inclusion of 
non-axisymmetric magnetic field terms. In particular, large numbers of distant 
negatively-charged dust grains, stable in a magnetic dipole, can be driven to 
escape by a more complex field. We trace the origin of the instability to over- 
lapping Lorentz resonances which are extremely powerful when the gravitational 
and electromagnetic forces on a dust grain are comparable. These resonances 
enable a dust grain to tap the spin energy of the planet to power its escape. We 
also explore the relatively minor influence of different launch speeds and the far 
more important effects of variable grain charge. Only the latter are capable of 
significantly affecting the micron-sized grains that dominate visible and infrared 
images of faint dust rings. Finally, we present full stability maps for Earth, 
Jupiter, Saturn, Uranus, and Neptune with magnetic fields modeled out to oc- 
tupole order. Not surprisingly, dust in the tortured magnetic fields of Uranus 
and Neptune show the greatest instability. 



1. Introduction 

When Voyager 1 encountered Jupiter in 1979, the discovery of the tenuous dusty ring 
system came as a complete surprise. Although the earlier Pioneer missions found some 
hints of a ring system, many thought that dust close to Jupiter would rapidly spiral in by 
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gas drag (lOwen et a/.lll979l ). Voyager 2 confirmed the existence of tfie ring, and found the 
tiny satelhtes Metis and Adrastea that orbit inside the classical Roche limit and are the 
most likely source of ring material. 

High speed impacts with small moons, as well as unseen large parent bodies, replenish 
the dusty rings with debris of all sizes. Similar sources of material for Saturn's tenuous inner 
D ring have not been found, though massi ve particles in narrow ringlets with e nhanced 



local densities could serve as these sources fIShowalter 



1996 



Hedman et al. 



20071 ). In both 



environments, dust ejected by impacts from parent bodies have essentially collisionless 
trajectories. As debris particles acquire electric charges through interactions with the 
plasma environment and solar radiation, the smallest reach significant charge-to-mass ratios 
and, as a consequence, experience strong electromagnetic (EM) forces as they orbit through 
the magnetic field of their host planet. 

For grains smaller than ~ l//m in radius, the EM force excee ds perturbations fro m large 



satellites, the planetary oblateness and solar radiation pressure (IHoranyi et al. 



19921 ). Even 



smaller dust grains may have orbits that are immediately unstable to either radial m otion 



if the grains are positively-charged ( iHamilton and Burns 



1993a 



vertical motion (both positively- and negatively- charged: 



Horanvi et al. 



iNorthrop and Hill 



1993ah . or 



19821 ). Many 



aut 
see 



lors have studied various aspects o f charged-particle dynamics. For a recent review. 



Jontof-Hutter and Hamilton 



(I2OI2I ). who derived analytic stability boundaries for the 
idealized case of grains with constant charge, launched at the Kepler speed in an aligned 
dipolar planetary magnetic field. As in that study, the boundaries between stable and 
unstable orbits are of particular interest to us here; these depend on the launch distance 
from the planet and the charge-to-mass ratio of an individual dust grain. 

The aim of this paper is to explore the sensitivity of th ese stability boundaries to m ore 



realistic situations. We relax the ideahzed assumptions of 



Jontof-Hutter and Hamilton 
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( 120121 ) above by considering i) non-zero ejecta speeds from the parent body, ii) higher-order 



magnetic field components, and iii) variable electric potentials on dust grains. We use 



well-s 


tudied dusty ring system ( 


Burns et al. 


1999; 


de Pater e 


al. 


1999; 


Ockert-Bell et al. 


1999; 


Brooks et al. 


2004 




Throop et al. 


2004 


: Showalter et al. 


200 


8; 


Kriiger et al. 


2009 


)■ 



After a detailed study of Jupiter, we then present stability maps for motion in the complex 
magnetic fields of Earth, Saturn, Uranus and Neptune. We begin by recapping stabilit y 



results for a simple dipolar planetary magnetic field from 



Jontof-Hutter and Hamilton 



(|2012h . 



2. Motion in an Aligned Dipolar Magnetic Field 

The charge-to-mass ratio for Kepler-launched grains can be conveniently described by 
the ratio of the force induced by the corotational electric field of the planet with gravity, 
given by 

= (1) 



( iHamilton 



1993a 



HAL 



Jontof-Hutter and Hamilton 



20121 ). Here, q and m are the electric 



charge and mass of a dust grain, qiq is the dipolar magnetic field strength at the equator, 
Rp and Mp are the planetary radius and mass respectively, VLp is the spin frequency of the 
planet, and G and c are the gravitational constant and speed of light. As a dimensionless 
independent variable, L^: accounts for all the relevant planetary parameters and avoids 
undue focus on the grain's size, shape, density and electric potential which are all poorly 
constrained. The sign of depends on the product qgiQ, and its value can easily be 
converted to a grain radius for specified grain properties. 



For large grains dominated by gravity and orbiting with semi-major axis a (the Kepler 
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limit), azimuthal, radial, and vertical motions have the same frequency 



rir 



(2) 



but for higher charge-to-mass ratios, these frequencies differ. As the charge-t o-mass ratio 



1993al ). and for an 



is raised, these frequencies slowly diverge from one another ( iHamiltonl 
aligned di polar magnetic field, explicit expr essions valid for all charge-to-mass ratios are 



available (IJontof-Hutter and Hamilton 



20121 ). As these expressions will prove useful for our 



current study, we reproduce them here. 

General motions in this problem can be conveniently separated into epicyclic motion 
about a guiding center which in turn circles the planet at an azimuthal angular speed Uc- 
Where radial epicycles are small on the scale of the grain's orbit, balancing the centrifugal 
force, the EM force and gravity yields an expression for Uc- 



2 , GMpL^ 
= uj^ + — 



(INorthrop and Hill 



1982 



Mitchell et al 



2003 



1 _ ^ 



GM„ 



(3) 



Jontof-Hutter and Hamilton 



2OI2I). Here and 



throughout, the subscript c refers to the guiding center of motion. The distance to the 
guiding center of motion, r^, is just the semi-major axis a in the Kepler limit. Note that for 
gravity-dominated grains (L* — )■ 0), we have u"^ = GMp/r^ = in agreement with Eq. [2l 
In the strong EM limit {\L^\ — ?■ ±00), Uc — > and the grains are nearly locked to the 
magnetic field lines. 



The radial or epicyclic frequency Kc satisfies 



gey 



(4) 



(IMendis et al. 



1982 



Mitchell et al. 



2003 



Jontof-Hutter and Hamilton 



20121 ) at the 



guiding center, where Qgc = qB/mc = n'^L^./Qp is the frequency of gyromotion. In the 
EM-dominated Lorentz regime, Kc = ^Igc- In the gravity-dominated Kepler regime, flgc — > 
and Kc — 7- Tic, the Kepler orbital frequency, as expected. 
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Most grains are radially confined, suffering excursions of 



' nl^-Q,Lm, + nL)+nl 
where r^, the gyroradius, is much smaller than ri, the launch distance. Here til = \/ GMp/' 



(5) 



and QgL = n^ L^ IVLp are the Kepler 



launch distance ( S chaffer and Burns 



requency and the gyrofrequency as dete rmined at the 



1994 



Jontof-Hutter and Hamilton 



20121 ). The epicyclic 



model fails only for posi tively-charged grains with L 



becomes very small (see 



Jontof-Hutter and Hamilton 



1, where the denominator in Eq. O 



2012 



, Fig. 6b). 



Finally, the vertical motion of grains with stable epicycles in the equatorial plane has 
frequency fif,, where 

3 



Pc 



2^ 



(6) 



( jJontof-Hutter and Hamilton 



Vtp is the azimuthal motion of the 



20121 1. Here, 0^ = 

guiding center in the frame rotating with the planet. Equation [6] is valid as long as 
'^g/^L « 1 which holds in both the Kepler (L* — t- 0) and Lorentz (L^, — )■ ±oo) limits. In 
the Kepler limit, all three terms are negligible and Vt^ — )■ ric, while in the Lorentz limit. 



only the last two terms can be ignored and Vt^ y^p ~ ^VLpUc + |n 



Where Vlb tends to zero, grains in the equatorial plane become locally vertically 
unstable. Equation |6] provides good agreement with numerical data on the location and 
charge-to-mass ratio of boundaries between vertically stable and unstable grains, with two 
important caveats. 

Firstly, in applying the epicyclic approximation, Eq. [6] assumes that radial motions are 
very small on the scale of the orbit (r^ << tl). In addition, Eq. [6] is averaged over one 
gyrocycle, so the epicyclic motion must occur on a much shorter timescale than any stable 
vertical oscillations >> ^b)- Both of these assumptions are easily met in the Lorentz 
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limit, but both lose accuracy as decreases, particularly for the positively-charged grains 
which become radially unstable as — ^ 1. 

Secondly, setting Qb = determines local, as opposed to global, vertical stability in 
the equatorial plane of the spinning planet and its aligned dipolar magnetic field. Local 
instability is a necessary condition for global instability (whereby grains collide with the 
planet at high latitude), but it is not always suf ficient. High- latitude restoring forc es often 



lead to stable high latitude oscillations (HLOs) fjJontof-Hutter and Hamilton 



20121 ). This 



class of orbits is more important for slow rotators like the Earth than it is at Jupiter or 
Saturn, but they do occur for the smallest grains inside 1.5Rp at Jupiter (Fig. [T]). 



With these two caveats in mind, we include the local and glo 



found from numerical integrations by 



3 al sta bility boundaries 



Jontof-Hutter and Hamilton 



(12012) for Jupiter with 



an aligned dipolar magnetic field model. Fig. [T] highlights these regions for a range of 
charge-to-mass ratios spanning four orders of magnitude and a suite of launch distances 
from Jupiter's surface to beyond its synchronous orbital dista nce, R.^yn, with grains all 



l aunch ed at the local circular speed of the large parent bodies. 



Jontof-Hutter and Hamilton 



(120121 ) also derived analytic approximations to most of the boundaries in Fig. [T] (see their 
Fig. 9). 



3. Jupiter 



We focus most of our attention on Jupiter as its magnetic fi eld has been wel. 



and is 


mown 


Dessler 


1983) 



mown out to octupole order (we adopt the 04 model of 



Acuna and Ness 



studi ed, 



1976 



Gauss, Qii = —0.664 Gauss, and hu = 0.264 Gauss; these can be combined to determine the 
dipole tilt angle: arctan ^ a/ {gfi + ) / Qioj = 9.6°. The (720 = -0.203 Gauss component can 
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Fig. 1. — Stability boundaries for a) negative (g < 0) and b) positive (g > 0) grains in 
an aligned dipole magnetic fi eld fo r Jupiter. This figure summarizes numerical data from 



Jontof-Hutter and Hamilton! ( 120121 ). with grain radii marked along the top axis correspond- 
ing to a —5V or +5V potential on a spherical grain of material density 1 g cm~^. The 
vertically-unstable grains depart from the equatorial plane and climb to high latitudes im- 
mediately after launch, ultimately colliding with the planet. Directly below this unstable 
region are grain trajectories with high latitude oscillations (HLOs). The radially unstable 
grains (g > 0) escape if launched outside synchronous orbit {Rgyn), or hit the planet if 
launched from within Rgyn- Two regions of high radial and latitudinal oscillations (HRLOs) 
abut the radial instability region. Within Rgyn, large inward radial excursions lead to vertical 
oscillations (roughly along magnetic field lines) which increase in amplitude until the grains 
strike the planet at high latitude. Outside Rgyn, near L^, = i, some grains experience HRLOs 
indefinitely. Finally, for g < 0, a curve traces grains that experience HRLOs following the 
2:1 resonance between the epicyclic (kc) and vertical (f2f,) frequencies. 
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be interpreted as a southward vertical offset to the dipole field. Four additional quadrupolar 
and seven octupolar terms are known, and the upcoming Juno mission will measure still 
higher-order magnetic field coefficients for the first time. In this section, we add various 
effects to a simple aligned dipole model to elucidate their importance. We begin with 
non-zero launch speeds in the frame of the parent particle, as likely occurs with impact 
ejecta. 



3.1. Varied Launch Speed 

Typical ejecta velocities from an impact are tens to hundreds of meters per second 
i n the rest frame of the pare nt body, in a cone centered on the impact velocity vector 



(Ide Pater and Lissauer 



20101) ■ Do these non-circular launch speeds significantly affect the 
stability of charged dust grains? To highlight the effect, we consider large initial velocities 
of 0.5 km s~^ in the prograde azimuthal (Fig. [2^) and radial (Fig. Hjo) directions. Even 
with these large speeds, we note that EM-dominated grains on the left side of the plots 
are hardly affected. The Kepler speed at = 2.0Rp is Vk =29.8 km s^^, while the local 
magnetic field lines rotate at QprL= 25.1 km s^^. The azimuthal impulse that we add, 
therefore, is only ~10% of the ^gc^g = 4.7 km s~^ gyrospeed and decreases the gyroradius 
Tg by a corresponding 10%. Although important, this effect is not noticeable on Fig. |2l 

On the other hand, grains in the Kepler regime experience large radial excursions 
following a launch impulse. For an azimuthal boost (Af^ > 0), we can solve for the 
semi-major axis a and eccentricity e from = a{l — e), and 

vl = GMp (---). (7) 

For Av^ << Vk, the radial motions extend outward from the launch position by 

2ae ~ AriAVfi^/vk ~ 0.134i?p for the parameters of Fig. [2^. Although the intermediate-sized 
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Fig. 2. — Radial range of motion for negatively-charged grains on initially circular orbits 
subject to a) an azimuthal (A^;^) and b) a radial (A^;^) velocity impulse. The small filled 
circles indicate particles launched on circular orbits with Av = 0, while the large open 
triangles denote those launched with Av — 0.5 km s~^. Synchronous orbit {Rsyn) is indicated 
by a dashed line; the other lines simply track the data points. 
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grains have the largest radial excursions in Fig. |2^, the grains in the Kepler limit are most 
strongly affected by a /S.v^ Icick. 

Figure Wp shows that a radial impulse produces a more modest radial range of motion 
than an azimuthal boost. In this case, the impulse is perpendicular to the velocity of the 
parent body. For the smallest grains, in the Lorentz limit, this has almost no effect on 
the motion perpendicular to the field lines, and is akin to altering the initial phase but 
not the size of the gyrocycle. As with the azimuthal kick, a radial impulse has the largest 
effect for the largest grains. To first order in Af^/ffc, the orbital energy and semimajor axis 
are unchanged. The range of motion is therefore centered on the launch distance and has 
magnitude 2ae = 2rLAvr/vk ~ 0.067 Rp for the parameters in Fig. [2]d. 

Figure [3] highhghts the effect of velocity impulses on the stability boundaries of Fig. [T] 
for positively-charged grains in three orthogonal directions: a prograde azimuthal impulse 
(Af0 = +0.5 km s""*^), a radial boost {Avr = +0.5 km s~^), and a vertical kick {Avz = +0.5 
km s~^). In each case, the orbital stability boundaries are only moderately affected by these 
changes; circular orbits are thus often a good approximation when considering stability. 
Only the azimuthal impulse appreciably affects the orbital energy, and hence shifts the 
radial stability boundary (larger grains on the right in Fig. |3^). In this case, the positive 

increases the Kepler orbital energy significantly, thereby preventing grains near the 
right-most radial stability boundary from falling into Jupiter. A negative Av^j, would 
destabilize grains near this boundary, permitting additional grains to fall to the planet. 
The left side boundary of the radially unstable zone is basically unaffected by all impulses, 
in agreement with Fig. [2l 

The vertical stability boundaries are moderately affected by the Af,^ (Fig. [3^) and 
Avz (Fig. [3b) initial impulses, but a radial impulse (Fig. [3b) has almost no discernable 
effect. Note that although in the Lorentz limit, the radial range of motion is too small to be 
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Fig. 3. — Stability of positive grains launched with a) an azimuthal launch speed Af^ = 0.5 
km s~^ faster than the local Kepler speed, b) Af^ = 0.5 km s~^, and c) Af^ = 0.5 km s~^, in- 
tegrated over 0.1 years. The solid curves are the numerically-determined stability boundaries 
of Fig. [1] where grains are launched at the local circular Kepler speed. The darkest regions 
denote grains that collide with the planet near the equatorial plane or escape. Moderate grey 
marks grains that collide with the planet at high latitudes (A > = 5°), while light grey 
shows HLO grains; locally unstable with vertical oscillations exceeding 5°, but bound glob- 
ally. White regions mark locally-stable trajectories. Here, radial stability (central regions) 
is only significantly affected by Af ^ while vertical stability (left-most regime) is affected by 
Ati^ and Af^. A radial impulse, Af^, has no noticeable effects. 
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significantly altered with a Av^ launch impulse, the change in the area of a gyroloop, which 
alters the magnetic mirror force, still noticeably affects the high stability boundary in 
Fig. [3^. Since At><^ > in Fig. [3^, the increased gyrospeed expands the gyroloop, leading to 
a stronger mirror force and hence a reduced region of vertical instability (Eq. |6]). Enhanced 
instability results for Af,^ < 0. A vertical impulse Af^ of either sign also leads to additional 
instability (Fig. For moderate values of in particular, the Af^ impulse causes the 
vertical instability region to dramatically expand near = 3, = 2i?p, and merge with 
the HRLO region of large radial and vertical oscillations (cf. Fig. [1]). Negatively- charged 
grains (not shown) are similarly affected by 0.5 km s~^ impulses. 

Overall, since the majority of real debris particles have much smaller speeds relative to 
their parent satellites than the 0.5 km s~^ considered here, we conclude that the stability 
boundaries are fairly insensitive to grain launch conditions. We note that new stability 
boundaries appropriate for non-circular initial orbits could be derived analytically using 



Hami 



2003 



tonian methods (INorthrop and Hill 



Jontof-Hutter and Hamilton 



1982 



Schaffer and Burns 



1994 



Mitchell et al 



20121 ). but as the effect is unimportant for our purposes. 



we turn instead to more complicated magnetic field geometries. 



3.2. Vertically Offset Dipole 

In this section we isolate the effect on orbital stability of Jupiter's dipole offset, 
modelled by the Qiq and g2o magnetic field terms. The maps in Fig. H] show that the offset 
field exacerbates the vertical instability for both negative and positive grains but has little 
effect on the radial stability boundaries. For the positive grains, the vertically unstable and 
HRLO zones overlap as in Fig. 

In the equator plane, the (720 niagnetic field is radial, and the corresponding v x B force 
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Fig. 4. — Stability of Kepler- launched negative grains (a) and positive grains (b) in a ver- 
tically offset dipole field for Jupiter, modelled with just the gio and (720 magnetic field 
terms, and integrated for 0.1 years. The curves indicate the numerical stability boundaries 
for the centered and aligned dipole configurations from Fig.[Tl In this map, the dark 
radially unstable grains that either hit the planet or escape at low latitude (|A| < = 5°). 
The moderately-grey regions are vertically unstable grains that collide with the planet at 
high latitude (|A| > A^), and the light grey regions show HLO stable grains that exceed Am 
in latitude, the same criteria that we have adopted for the aligned dipole magnetic field of 
Fig. El 



- 15 - 



is vertical. Thus the (720 magnetic field term primarily adds an additional vertical force, 
thereby expanding the vertical instability. In fact, the stability map in Fig. [Do resembles 
Fig. [3t, which modelled a vertical impulse on the grains at launch; comparison of the two 
figures shows that the inclusion of g2o is a far more important effect. Note also that, unlike 
the effect of a Avz impulse, the offset dipole is effective at destabilizing grains near the 
planet, causing a significant vertical bounce oscillation and significantly expanding the 
region of global vertical instability. 

3.3. Tilted Dipole 

In testing the effect of a tilted dipole field, we include the 5^10, gn and hu magnetic 
field terms in our numerical models, setting g2o and all higher order terms to zero. Since the 
magnetic and gravitational equators do not coincide, we consider two separate equatorial 
launch phases: (i) (pQ = 0°, the ascending node of the magnetic equator on the geographic 
equator and (ii) 0o = 90°, where the magnetic equator reaches its highest northern latitude 
of 9.6°. At this launch phase, many grains can reach latitudes ~ 20° north and south of the 
equator, even if their trajectories are stable. Our stability results for negatively-charged 
grains are plotted in Fig. O 

Although slight differences with launch phase are apparent. Figs. [5^ and|5]D are quite 
similar. Jupiter's tilt is a stronger effect than its offset (see Fig. H^), extending the vertical 
instability boundary significantly outwards and close to Rsyn- The dramatic outward 
expansion of the vertical instability can be understood as follows. For an aligned dipole, as 
synchronous orbit is approached, both the velocity relative to the magnetic field and the 
electromagnetic forces tend toward zero. Furthermore, as the velocity is azimuthal and the 
field is vertical, the direction of the weak EM force is entirely radial. For a tilted dipole, 
however, the magnetic field lines cross the equator plane with a radial component, causing 
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Fig. 5. — Stability of negative grains integrated over 0.1 years in a tilted dipole field 
for Jupiter, with the launch longitude at two locations; a) where the magnetic equator 
crosses the planetary equator plane (0o = 0°), and b) where the magnetic equator reaches its 
highest geographic latitude (0o = 90°). As before, the curves are the numerically-determined 
stability boundaries for the aligned dipole case from Fig. [1^. The grey scale is similar to 
Figs. 13] and S) dark points are grains that collide with the planet at low latitudes, the 
moderately-grey region denotes grains that were vertically unstable to collide with the planet 
at latitudes higher than = 20°, the lightest grey marks trajectories that were excited to 
higher latitudes but remained bound, and the white areas represent grains that were locally 
stable. The only difference from Fig. [3] and H] is that here we define low latitude to be 
lA^I < 20° rather than 5°. 
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a substantial v x B vertical force as with the offset dipole. These forces push particles out 
of the plane along field lines, leading to an expansion of the instability zone nearly to Rgyn- 
Interestingly, the inner boundary at {L^ = —50, r^/Rp = 1.4) is far less affected. 

One key difference in the two panels of Fig. [5] occurs for high L^, along ri/Rp = 1; 
launching at the node (0o = 0°) leads to collisions while launching at 0o = 90° does not. 
This difference is due to the curvature of the field lines in a dipole. In an aligned dipole 
magnetic field, stable mirror motion causes EM-dominated grains to oscillate about the 
magnetic equator, whereby the turning points or mirror points confine the latitudinal range 
of the grain. Launching from 0o = 90° in the tilted magnetic field ensures that the launch 
point is at one of the mirror points, and this vertical turning point is relatively close to 
Jupiter. Thus grains launched near IRp in Fig. [5]d initially move radially outward and do 
not collide with the planet. By contrast, for grains launched at the node where (po = 0°, the 
mirror points are necessarily closer to the planet than the launch distance and, accordingly, 
we see that grains launched within l.OGRp are forced to collide with Jupiter in Fig. [5^. 

Another, more subtle difference between Figs. andOD, is that, grains launched at 
the node (0o = 0°) are slightly more stable close to Rgyn than those launched at 0o = 90° 
e.g. at {L^ = —100, r^/Rp = 2.21). In fact, both the inner and outer vertical stability 
boundaries are shifted slightly outwards for 0o = 90° compared to launches at 0o = 0°. 

Positively-charged dust grains are similarly affected by the addition of the dipole tilt. 
As with negative grains, we present two launch phases in Fig. El and find differences in 
orbital stability similar to those already discussed for Fig. |3 In particular, the azimuthal 
dependencies for highly-charged dust grains >> 1) are almost identical for both 
negative and positive charges (Fig. |6]). As in Fig. [5^, grains launched immediately above 
the planet are unstable for 0o = 0° (Fig. 

For all launch longitudes, the vertical instability expands greatly outwards, nearly 
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Fig. 6. — Stability of positive grains integrated for 0.1 years in a tilted dipole field for 
Jupiter, with two different launch longitudes: a) 0o = 0°, where the magnetic equator crosses 
the planetary equator plane, and b) (pQ = 90°, where the magnetic equator reaches its highest 
geographic latitude. The solid black curves mark the stability boundaries for the aligned 
dipole case from Fig. [TJo. As in Fig. |5l the darkest points denote grains that strike the planet 
at low latitude (|Am| < 20°, within twice the tilt angle), the moderate grey marks grains 
that strike the planet at high latitudes, the lightest grey marks grains that remain bound 
between high-latitude mirror points, and the white area represents grains that are vertically 
stable. 
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to synchronous orbit but very little towards the planet. As in Fig. O both the inner 
and outer vertical stability boundaries in the Lorentz limit are shifted slightly outwards 
for the 00 = 90° launches of Fig. [Hb compared to the 0o = 0° launches of Fig. [B^. The 
two boundaries have slightly different explanations. For the outer boundary near Rsyn, 
launching at (po = 90° allows the initially larger vertical electromagnetic forces to drive the 



grain to higher latitu des where the gravity of the p 



and cause instability ( jJontof-Hutter and Hamilton 



anet can overwhelm the centrifugal force 



2012). Near the inner vertical boundary, 



however, grains launched at the node 00 = 0° have a higher latitudinal range and are 
shghtly less stable. 

The radial stability boundaries are largely unaffected by the tilt in the magnetic field, 
although some slight differences are evident to the left of the radial instability region. Note 
the subtle difference along the left-most radial stability boundaries between Fig. and 
Fig. [6h, where grains launched at (po = 90° inside Rgyn, are slightly more stable than those 
with 00 = 0°. Outside synchronous orbit, however, the reverse holds true. This is most 
easily understood as an overall outward shift of the instability region from 0o = 90° to 
00 = 0°. Thus 00 = 90° grains behave almost exactly like 0o = 0° grains that have been 
launched a bit further from the planet. Note that this difference with launch phase was also 
seen for the negative grains with tl ^ Rp (Fig. E]), and the explanation is the same. 

Until now we have only considered instabilities that remove a grain typically within a 
few hours. However, the tilted dipolar field causes further instabilities acting over weeks 
to months, and over a greater range of launch distances than the aligned dipolar case. We 
explore these longer-term effects below. 
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3.4. Resonant Effects in a Tilted Dipole Field 



In an aligned dipole field, it can be shown that negative grains outside Rsyn are 



permanent 



y confined between their launch di stance and Rsyn', they are energetically unable 



to escape (jJontof-Hutter and Hamilton 



2012I ). The tilted dipole field however, permits 



radial motion away from Rsvn., and actually enables some negative grains to depart the 



planet, as was first seen by 



HamiltonI fll996h . 



In Fig. [3, we show the stability maps for Jupiter modelled with the gio and gn 
magnetic field components for dust grains with both negative and positive charges whose 
trajectories were integrated for one Earth year. We seek to highlight the motion away from 
synchronous orbit for the negative grains, and towards Rgyn for the positive grains, motions 
precluded by a simple aligned dipolar magnetic field. 

Within synchronous orbit, the negatively-charged grains of Fig. [7^ shows the same 
short-term instabilities seen in Fig. [Sb- Notice in Fig. [7^, however, the large fingers of 
instability outside synchronous orbit. These features trace grains that suffer significant 
and unusual motions away from Rgyn, and the largest one points tow ards rr./Rp = 3.55 



which happens to be the location of the outer 1:2 Lorentz resonance (ISchaffer and Burns 



1987 



1992 



Hamilton 



1994 



Showalter et al. 



20081 ). The other fingers point towards other 



Lorentz resonances. Importantly, the highly- detailed dark grey structures within these 
fingers (Fig. [7^) indicate significant numbers of negative grains that actually escape from 
Jupiter within one year. 

These results are important for the escape of du st from the lo plasma torus, the rn ost 



likely source of the jovian high-speed dust streams (IHoranyi et al. 



1993b 



Graps et al. 



2000l ). Dust streams are comprised of radially-accelerated positively-charged dust grains. 



although in the plasma torus itself, dust-grain elec tric potentials are likely to be negative 



even in sunlight (iBagena. 



1994 



Kriiger et al. 



20031 ). Lorentz resonances can provide a rapid 
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Fig. 7. — Destabilising effects of a tilted dipole magnetic field for Jupiter, with launch at 

00 = 90°. Integrations are for one Earth year, and we show a greater radial range than in 
Figs. |5] and El Note that the greyscale that we use here is also different from the earlier 
figures. Here, the darkest region denotes grains that either escape or crash into Jupiter within 

1 year of launch. For the negative grains, the light grey marks grains with radial motions 
away from Rsyn (in the direction opposite that expected for gyromotion), by at least OMri, 
revealing the destabilizing effect of the tilted magnetic field. For positive grains, the light 
grey indicates trajectories with radial motions towards Rsyn of at least 0.02ri:,. As always, 
white indicates stability. 
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escape mechanism for negatively-charged grains launched in the plasma torus. Once dust 
grains are free of the torus, charging currents become positive and the grains are accelerated 
outwards to escape. Our modelling in Fig. [7^ actually understates the importance of this 
mechanism, as we force the potential to remain negative far outside the actual boundaries 
of the plasma torus. 

For the positive grains. Fig. [7)d shows an increased number of grains that are unstable, 
compared to Fig. [6]d. The most striking difference is that in Fig. [Tb, there are rough 
patches of additional radial instability just outside Rgyn near {L^ ^ 0.2, ri/Rp = 3). These 
unstable patches transition smoothly to become thin tracks of bound grains with excited 
radial ranges in the Kepler regime which, like the negative grains in Fig. [T^i, point towards 
Lorentz resonances which occur for discrete integer ratios of the planetary spin and Kepler 
orbital frequencies. 



Accordingly, we look to extend t he concept o 



the Kepler limit by < 


authors 


Hamilton and Burns 


1993b: 



Hamilton 



Burn 



s et al. 



Lorentz resonances, (much studied i n 



1985 



Schaffer and Burns 



1987 



1992 



19941 ). to cover the entire range of charge-to- mass 
ratios. To determine the location of these Lorentz resonances in general, we start with the 
resonant equation 



^ = Auc + BVLr, + CVL^ 



node 



peril 



(8) 



199J). 



where the coefficients A, B, C and D are integers that must sum to zero (iHamiltonl 
Here ujc is the orbital frequency of the guiding center, Vtp is the planetary spin rate, O^node is 
the precession rate of the ascending node, Wperi is the precession rate of the pericenter, and 
^ is the resonant argument; ^ equals zero at a Lorentz resonance. 

Equation [8] is completely general and valid for all charge-to- mass ratios if we are careful 
to rewrite zuperi and Vtnode in terms of our fully general frequencies from Eqs. [3] - El The 
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precession rates are simply differences between fundamental frequencies: 
and 

^node = <^c~ (10) 

where Kc, the epicyclic frequency of motion, is negative by convention for retrograde 
epicycles. Recall that in the gravity limit, Uc = l^d = = and hence Wperi = ^node = 0, 
as expected. As an illustration, we focus on radial resonances for which C = 0. 

Setting Eq. [8] to zero and using Eq. [9] to eliminate Wperi, we find: 

- Bujc + BQp- D\kc\ = (11) 

or 

B4>, + D\n,\=0, (12) 

since in the frame corotating with the magnetic field, the azimuthal frequency of the guiding 
center is given by (pc = — ^p- 

Equation [T2] shows that a Lorentz resonance affecting radial oscillations reduces to 
a simple ratio between the epicyclic frequency \kc\, and the motion of the guiding center 
relative to the rotating magnetic field {(pc)- Our approach thus shows how to extend 
classical Lorentz resonances to remain valid at arbitrary charge-to-mass ratios. 

In Table 1, we show select Lorentz r esonances f or all charge-to-mass ratios and their 



driving magnetic field terms, taken from 



HamiltonI ( 1l994j ). In the Kepler limit, these 



Lorentz resonances act to slowly increase eccentricities and/or inclinations, destabilizing 
trajectories over many orbits. The resonances that include multiple instances of zuperi or 
flnode, such as the three 1:3 and 1:4 resonances in Table 1, are weaker since their strengths 
in the Kepler regime depend on higher powers of the small quantities e (eccentricity) and 
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5-field terms 


Res. name 


ri {L, ^ 0) 


Res. frequency 


Corotating form 


9n 


1:4 


5.64 Rp 


4c^c 2'ZZ7pg|.^ ^node 


0c + 2 + rife 


9ii 


1:3 


4.66 Rp 


3u^(. rip ZUpg^j-i ^node 


0c + «c| + ^^6 


9n 


1:2 


3.55 Rp 


2cUc ^node 


4>c + rife 


921 


1:3 


4.66 Rp 


3u^(. 2'rf7peri 


0C+ 2|/«c| 


921 


1:3 


4.66 Rp 


3u^c '^^node 


0c + 2Q6 


921 


1:2 


3.55 Rp 


2cJc ^peri 


0c + «c 


922 


2:4 


3.55 Rp 


^UJc 2 rip "^peri ^node 


20c + \Kc\ +^6 


922 


2:3 


2.93 i?p 


3u^c 2Qp rijjQ^g 


20c + % 


922 


2:1 


1.41 


2f2p + ^node 


20c - 



Table 1: Selected Lorentz resonances (column 2), driven by the magnetic field coefficient in 
column 1, appear, for small charge-to-mass ratios, at the locations given in column 3. The 
resonance frequency is given in its most general form (column 5) and in a second form most 
useful when gravity dominates (column 4). 
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i (inclination). At higher charge-to-mass ratios, however, all of these resonances increase 
in strength, and their effects on grain orbits occur on much shorter timescales than in the 
Kepler regime. Some negative grains at 6.0-Rp in Fig. [7] escape in as little as a few days. 

Figure [8] overlays the strictly radial Lorentz resonances of Table 1 on the stability map 
for a tilted dipolar field (the data from Fig. [7^). For negative grains outside synchronous 
orbit, the Lorentz resonances curve upwards directly into the region of escape for increasing 
L=K. This occurs because the epicyclic frequency \kc\ increases rapidly with (Eq. Hj); 0c 
must also increase to maintain a given resonance (Eq. [T2l) . Since 0c increases away from 
Rsyn, remaining in resonance as \kc\ increases necessitates a greater launch distance from 
synchronous orbit. Although these curves are determined from frequencies that are strictly 
valid only for an aligned dipole field, they nevertheless show an impressive match to our 
data, despite the more complex magnetic field. 

For the positive grains, the resonant tracks in Fig. [8]d begin at the same locations 
in the Kepler limit as for negative grains, but they curve towards synchronous orbit as 
L^, increases. On the right side of the short-term radial instability of Fig. [H all resonant 
solutions converge to a single point at (L^. = 2 — \/3, = Rgyn)- This is the point at 
synchronous orbit wh ere Ivl = 0c = 0, and grain orbits are locally unstable in even a simple 



aligned dipolar field ( iJontof-Hutter and Hamilton 



20121). The convergence to synchronous 



orbit as L^. increases occurs because for positive grains with << 1, \kc\ decreases as 
increases (Eq. Hj), and so 0c must decrease as well (Eq. [12]), driving distances towards Rgyn- 

The Lorentz resonances destabilize the motion of grains, and hint that a non-axisymetric 
field allows the negative grains to tap into planetary rotation, to make escape energetically 
favorable. The detailed structure in the stability map of Fig. |8l including the escaping 
negatively-charged grains, is due only to the effects of qiq and gu. The first-order theory, 
however, can only explain the three 1:N resonances (Table 1) and not the instability of the 
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Fig. 8. — The tilted dipole stability map of Fig. [7] with theoretical curves = for radial 
resonances only from Table 1 superimposed as solid curves. The theoretical curves fall atop 
the instability "fingers" seen in Fig. [TH, and the trails of points from Fig. Uh, attesting to 
the accuracy of the theory. The 2:1 resonance between epicyclic frequency \kc\ and vertical 
Qfy motions is also shown in the upper panel, as a dashed curve. In the lower panel, the open 
circles mark the points (L^: = 2 ± v^, ri/Rsyn = 1), where (pc = = and all resonant 
tracks for positive grains converge. For large L*, all radial resonances lie nearly atop one 
another in panel b). 
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2:3 and 2:1 resonances, which nevertheless are definitely present in Figs. [7] and El We will 
return to explain this discrepancy shortly. 

In addition to Lorentz resonances of the type shown in Eq. [HI there are also resonances 
between the dust grain's radial and vertical motions, analogous to the Kozai resonance 
experienced by highly-inclined orbits. The dominant resonance of this type satisfies: 
LOc — 2Qnode + ^peri = 0, such that \k,c\ = 2i7fc. This 2:1 resonance between radial and 
bounce motions is the strongest of its type since during one bounce period, north-south 
symmetry ensures that the dust gr ain experiences two cycles in magnetic field strength 



(jJontof-Hutter and Hamilton 



2012). The resonance track also passes close to the high 
charge-to-mass boundary of the resonant structure in Fig. [8^. We turn now to investigate 
the effects of the individual asymmetric quadrupolar magnetic field terms which should also 
power resonances (Table 1). 



3.5. Quadrupole Terms 

In this and the following sections, we focus on the escaping negative grains outside 
synchronous orbit, because these escapes are the most fundamental new effect added by a 
non-axisymmetric magnetic field. In the stability maps of Fig [9], we isolate the effects of 
(721 and (722 to highlight their respective Lorentz resonances, as compared to the stability 
boundaries of the tilted dipole from Fig. [7^ (solid curves). For stable negatively-charged 
grains in an aligned dipole field, radial motion is always confined between the launch 
distance and synchronous orbit. Thus, as in Fig. [7^, the light grey data in Fig. [9] trace 
where grain trajectories show motions away from Rgyn that are significant on the scale of 
the launch distance. 

The 5^21 and g22 terms studied in Figs. [9^ and [Qjo clearly cause less overall instability 
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Fig. 9. — Quadrupole-order terms isolated in stability maps for negatively- charged dust: 
To QiQ we add just (721 in panel a) and just (722 in panel b), with a 1-year simulation for 
each dust grain. The light grey marks stable grains whose radial excursions away from Rsyn 
exceeded 4% of the launch distance, while dark grey indicates collision just as in Fig. [7^ 
andlHk- The dark curves mark the envelope of instability when just the gio and gn coefficients 
are included, from the numerical data of Fig. [7^. 
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than Qii in Fig. [7^. Consider first, motion within synchronous orbit. The (721 term (Fig. [9^) 
is nearly as effective as gu in inducing vertical instability, and in fact is better able to clear 
out the region just above the planet's cloudtops. The effect of the (722 term (Fig. [Hb) is 
similar to, but typically weaker than, with one important exception. Note the long, 
horizontal feature extending towards the 2:1 inner Lorentz resonance in the Kepler regime 
(small L*). These stable grains are strongly stirred by the 2:1 vertical Lorentz resonance 
excited by (722 (Table 1). A similar effect can be seen on the right-hand side of Fig. [9^; a 
small trail of points near (r^/i^p ^ 1.05, L^, ^ —0.03) hints that q ^i cau ses a weak 2:1 inner 



HamiltonI (Il994j ) predicts that (732 



Lorentz resonance. However, the first-order theory of 
rather than (721 should excite this resonance! 

Outside Rsyn, the situation is more straightforward. Figure [9^ shows that the (721 term 
strongly excites the 1:3 and 1:2 Lorentz resonances, as expected from Table 1. Notice that 
the 1:2 resonance is significantly stronger than the 1:3 resonance because, in the Kepler 
limit, the former has a strength proportional to the small orbital inclination i while the 



latter's strength de pends on t 



orbital eccentricity ([Hamilton 



le pr oduct of two small quantities e and i, where e is the 



1994j ). Interestingly, in one way the (721 term has a stronger 



effect on radial motion than the gu term, extending the 1:2 and 1:3 re s onanc es further into 



the Kepler limit. This can be understood from Table 1 and 



Hamilton! (119941 ). which show 



that the (721 term should naturally excite both of these resonances. 

Figure [9]d shows two main features from the (722 term outside Rgyn, which tend towa rds 



mi), 



1:2 and 2:3 in the Kepler limit. The (722 term excites vertical motions (iHamiltonl 
which are not traced directly in Fig. |9l but which clearly couple to radial motions. This 
causes the outer 1:2 and 2:3 resonances seen in Fig. [9]d, as well the strong inner 2:1 
resonance that reaches far into the Kepler regime. A glance at Table 1 shows that (722 
excites a first-order 2:3 inclination resonance and a second-order mixed and therefore weaker 
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2:4 resonance, accounting for the differing responses of grains near these resonances visible 
in Fig. [Hb- 

We are left with a few paradoxes. First, how does the gu magnetic field term excite 
the 2:3 and 2:1 Lorentz resonance? And how does the (721 term excite the inner 2:1 
resonance? To answer t hese question s, we require a second-order expansion of the Gaussian 



perturbation equations (jPanby 



19881 ). for the electromagnetic force. Th e first- order Fourier 



Hamilton! (119941 ) for each 



series expansion in the small parameter was obtained by 
magnetic field coefficient by treating the orbital elements as constants. To extend this 
to second order, we take the Fourier series first-order solution for each orbital element 
and insert it on the right-hand side of the Gaussian perturbation equations. Simplifying 
requires identities for the product of two trigonometric functions and we end up with 
second-order Ll corrections to the time rates of ch ange of the orbit al elements. Thus 



Hamilton! (jl994j ) is augmented by a 



the power in each resonant frequency in Table II of 
second-order correction. Calculating the strength of these corrections is a straight-forward 
but unenlightening exercise which we do not undertake here, as the calculation is clearly 
invalid for L,, > 1 when the third- and higher-order terms cannot be ignored. Indeed, the 
very concept of orbital elements also breaks down for > 1 when electromagnetism is no 
longer a small perturbation to gravity. 

Instead, we explore the form of the corrections and show how magnetic field c oefficients 



can excite resonant terms other than those shown in our Table 1 and in Table II of 



Hamilton 



(I1994J ). Consider first the gioQu simulation of Figs. [7]and|Hl To second order in L^,, this 
combination of coefficients excites two relevant 2:3 resonances: \1/ = — 2Qp — Wp^ri and 
\1/ = ?>ujc — 2i7p + vuperi — 2Qnode, botli with amplitude proportional to Ll{Rp/r)^gl^ei'^. 
In addition, a 2:1 resonance, \l/ = — 2f2p + Wperi is also excited, with amplitude also 
proportional to Ll{Rp/r)^gf^ef . These resonances show up in Figs. [7] and [8] near the planet 
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where [Rp/r) is relatively large, near instability boundaries where e and/or i are large, 
and for particles where itself is relatively large. Comparison of the data in Fig. [7)d 
with the corresponding curves in Fig. IHb shows that these resonances are weaker than the 
already-discussed first-order resonances, as expected. 

In a similar manner, the second-order theory shows that g2i also drives the 2:1 
resonance with frequency \E' = Wc — 2Qp + Wperi and amplitude Llgli^Rp / r)^ (Fig. |9K). 
he r est of Figs. |9^ and [9]d appear to be well explained by the linear theory of 



Hamilton 



( 1l994j ). which predicts both pairs of instability outside i?, 



syn- 



The second-order corrections, however, should excite the resonant frequency 
\I/ = — 4i7p + Wp^ri — '^^node with amplitude proportional to Ll{Rp/ r)^ g22eP ; the tiny 
weak feature just below the 2:3 track and near the center of Fig. ^ may be due to this 
resonance. Furthermore, the frequency \E' = Su^ — 2Qp — Wperi with amplitude proportional 
to Ll{Rp/r)^g2ie^, should excite particle motions in Fig. [9^, but no evidence for these 
motions is seen. This may be due to the fact that horizontal re sonances driven by g2i are 



intrinsically weaker than the vertical resonances driven by g22 ( iHamilton 



19M). In any 



case, given the strong drop in the strength of second-order corrections with distance, their 



effects outside Rsyn are minimal. 



3.6. Realistic Full Magnetic Field Models 

Figure [TOk combines the effects of all dipolar and quadrupole terms for negative grains. 
Within synchronous orbit, all grains in the Lorentz limit are now unstable, as are all grains 
within the gu envelope. Furthermore, (722 powers a 2:1 vertical resonance with frequency 
\I' = Uc — 2Qp + Clnode, seen as the dominant horizontal feature extending from the instability 
region (inside Rgyn in Fig. [TOb). Adding the octupole term 5(32 strengthens this feature by 
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exciting a 2:1 radial resonance \1/ = — 2^2^ + vjperi (Fig. [TOb). These resonances have the 
most dramatic effect on large dust grains near the planet, due to rapid decline in magnetic 
field strength with distance for higher-order terms. The (733 octupole term adds a spike of 
instability at the 3:2 resonance as well (\E' = 2uJc — 3Qp + ^Inode)- 

Outside Rsyn, a large region of escaping negative grains exceeds the sum of the effects 
of gii (Fig. [7^), (721 (Fig. [9^) and (722 (Fig. [9]d), although the main resonant tracks are easily 
identified. In particular, a huge swath of grains centered on = 5 at r^/i^p = 5 escapes 
here, but is bound for the simpler field geometries of Figs. Uh-, [9k, and [9b. 

Adding the octupole magnetic field coefficients (Fig. [TOb) presents only subtle 
differences from the quadrupole model of Fig. [TOk outside Rsyn- In particular, the locations 
of the resonant tracks appear to be unchanged. The three narrow fingers in the center of 
Fig. [TOb . however, are noticeably more prominent than the corresponding structures in 
Fig. [TOk . The outer 2:3 resonance is driven by the g22 term (\E' = Suc — 2Qp — Qnode), but 
the 3:4 and 4:5 resonances cannot be excited by quadrupole terms in the linear theory. The 
3:4 resonance is driven by the (733 coefficient, but also by a second-order term proportional 
to gug22- Both are active in the lower plot, while only the latter affects the upper plot. 
Similarly, the 4:5 resonance is excited by the non-linear (^22 term (both plots) and by the 
(711(733 term (bottom plot only). As always, when multiple resonances are active, chaos 
ensues and escape becomes more likely. Note that these differences between Figs. [TOk 
and [TOb are confined within ~ 4i?j,, due to the rapid radial weakening of the high-order 
magnetic field terms. 

In general, we see numerically and analytically that Lorentz resonances widen in 
strength as increases. This causes the resonances to overlap and destabilize most of the 
grains near = — 1 if grains are launched beyond the immediate vicinity of synchronous 
orbit. As |L*| increases, higher-order dependencies on the charge-to- mass ratio permit even 



- 33 - 




100 10 1 0.1 0.01 

-L* 



Fig. 10. — Resonant quadrupolar and octupolar effects in stability maps for negative 
dust grains: To the tilted dipole model of Fig. [7^, shown here as the solid curves, we add 
all quadrupole terms (panel a), and all terms out to octupole order (panel b). The darkest 
regions mark grains that collide with the planet or escape during a 1-year integration. The 
lighter grey indicates stable grains whose radial excursions away from Rsyn exceeded 4% of 
the launch distance, as in Figs. [7^, [8^1 and [91 
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more resonances to emerge and vie for control of dust grain dynamics. 

In Fig. [m the vertical and radial Lorentz resonances for negative and positive grains 
are atop the stability map for Jupiter's full magnetic field modeled out to octupole order. 
For the negative grains of Fig. [Tib , as increases going from right to left, all the radial 
resonances diverge rapidly from Rsyn- Most of the vertical resonances however, diverge 
from synchronous orbit more slowly as \L^,\ increases, and in the Lorentz regime these pile 
up on the vertical stability boundary inside Rsyn (Fig. [T^), where fif, — and hence, by 
the resonant condition, <pc 0. The combined effects of many vertical resonances near this 
boundary destabilizes all grains in the Lorentz regime out to synchronous orbit in Fig. [TTh . 

For the positive grains of Fig. [TTb . all radial resonances converge on the two locally 
unstable points along Rsyn- In the Lorentz regime, the curve outside synchronous orbit 
satisfies Kc — > 0. This is further to the left of the stability boundary for an aligned dipole 
(Fig. [TJj), ensuring that resonances pile up and further destabilize grains with the additional 
magnetic field terms in the Lorentz regime. Physically, it means that smaller grains are 
more likely to be expelled for a particular positive electric potential than calculated using 
the aligned dipolar approximation. For the larger grains, in the Kepler regime, the vertical 
Lorentz resonances asymptote near the L^, = 1/2 boundary where the guiding center 
distance rapidly increases and the bounce frequency — )■ 0. 

For both positive and negative grains where \L^\ << 1, the outer 1:3 radial and vertical 
resonances coincide i mplying zjppri = ^ nnd,-.- This is indeed the case to first order in L^: as 



HamiltonI (jl993al ). The result can also be obtained from our Eqs. [31 [H 



was first deduced by 
and [61 

In the Kepler regime, N:N+1 resonances pile up at synchronous orbit (like the 5:6 
resonance marked in Fig. [TTb): for higher N, these are driven by g^^- and Qnn-i magnetic 
field terms beyond the octupole model that we have considered here. Thus our numerical 
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Fig. 11. — Theoretical resonance curves over all cliarge-to-mass ratios superimposed on a 
stability map of Jupiter's full magnetic field for a) negative grains (from Fig. [TOb) and b) 
positive grains. The dark grey marks dust grains that escaped or crashed into Jupiter during 
the 1-year integration. The faint grey points denote grains that experience radial motions 
away from synchronous that exceeded 4% of the launch distance for the negative grains, 
and towards Rgyn by at least 2% of for the positive grains. The thick bold curves mark 
radial Lorentz resonances, and the thin curves track the vertical Lorentz resonances. At 
synchronous orbit (t>r = 0, a nd th e white circles mark the local stability threshold from 



Jontof-Hutter and Hamilton! (120121 ) where Kc 0. 
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model of Jupiter's magnetic field is incomplete and the inclusion of higher-order terms 
would lead to some additional escapes. Given the strong radial dependence of higher-order 
magnetic field components, we expect changes to be limited to regions close to the planet 
and at high L*, just like the differences between Figs. [TOb andfTOb. Nevertheless, we eagerly 
await the improved magnetic field model that the Juno spacecraft will soon provide. In the 
meantime, we now relax our assumption of a constant grain-charge. 



4. Variable grain charge in an aligned dipolar field 



The electric charging of a circumplanetary dust grain is a function of the plasma 
environment, the flux of solar radiation, and the physical properties of the grain itself. 
Since the nature of the grains and their plasma environment are poorly constrained, the 
motion of any particular grain with a varying charge is highly model dependent. Our goal 
in this section is not to pick the best model for a given situation, but rather to elucidate 
the physics of orbital changes driven by charge variations. Possibly the simplest non-trivial 
model which nevertheless, must occur in circumplanetary ap plications is the shutoff o f 



the photoelectric current during planetary shadow passages (iHoranyi and Burns 



199ll) 



This effect will be present even if all other model-dependent charging effects are absent. 
The azimuthal asymmetries that the shadow induces has a profound effect on dust grain 
motions, as we shall soon see. 

Returning to our aligned dipolar magnetic field model, and assuming that the plasma 
distribution is perfectly axisymmetric, the effect of the planetary shadow transit is to 
introduce a strong azimuthal asymmetry in the charging environment for a dust grain. In 
the shadow, the photoelectric effect of sunl ight is absent, and intera ctions with the plasma 



cause a net negative charge on the grains fISchaffer and Burns 



19871 ). In the sunlight, by 



contrast, equilibrium typically favors a slight positive electric potential. Fig. [T2] highlights 
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Fig. 12. — The charge response of a dust grain in plasma and sunhght depends on grain 
size. Here we show a) 10 fim and b) 0.0527 fim grains, launched at 1.4 Rp at local noon, 
with an aligned dipole magnetic field for Jupiter, in a uniform plasma with density Ue = 1.4 
cm~^ and temperature Tg = 10 eV. The grains initially carry no charge. Each panel shows 
the grain's radial trajectory and instantaneous electric potential $ in Volts. The large grain 
experiences three 1-hour long shadow passages during which the charge decreases, while the 
smaller grain has a single 3-hour eclipse, and much more sluggish changes to its electric 
potential. 
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the differing effect of charge variation with grain size, given identical launch distances. 
Here, we have chosen two example grains launched inside synchronous orbit at a location 
that avoids the short term vertical instability of Fig. [T] 

Firstly, we note that the equilibrium charge on a sunlit dust grain is 2.75 Volts, whereas 
in the shadow, the equilibrium electric potential is -27 Volts. The large grains reach their 

ndeed, the charge respon se 



20101 ) 



equilibrium potential far more rapidly than the smaller grains. 

is typically inversely proportional to the size of the dust grain (iHoranyi and Juhasd 
For the larger grain in Fig. [T2k . the increasing amplitude of the radial oscillations is caused 
by the fact that charge variation repeats each dust grain orbit, thereby resonating with 



the epicycl ic frequency for grains in the Kepler regime. T 



resonance jHoranyi and Burns 



1991 



Hamilton and Kriiger 



lis is the destabilizing shadow 



20081 ) which we will find strongly 



affects our stability map. 

In Fig. [T2b . the smaller dust grain does not have enough time to reach charge 
equilibrium during its three-hour shadow passage. This dust grain experiences stochastic 
kicks both radially inwards and outwards from its launch distance; the grain eventually 
becomes vertically unstable and crashes into the planet at high latitude after just 14 hours. 
Each kick in the guiding center distance occurs when the electric potential on the dust 
grain is ~1 Volt, when the instantaneous value L^, places the grain near the left-most radial 
stability boundary of Fig. [T]d. When the potential is higher than 2 Volts or negative, the 
grain experiences stable radial oscillations. The decreasing amplitude of these oscillations 
with time is due to the grain reaching higher charge-to-mass ratios and hence 

experiencing tighter gyrations. After several random steps in due to the grain's periodic 
encounters with the radial instability, it moves into the vertical instability zone (Fig. [1]) and 
is lost to Jupiter. 



In Fig. [13] we present stability maps for a large range of grains sizes from O.OOl/im 
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to lOyum, over a broad range of launch distances, to test the effect of charge variations 
on orbital stability. In these models, dust grains are free to vary their charge as the 
environment allows, both with the effects of the shadow present (Fig. [T^). and explicitly 
ignored (Fig. [T3b). 

Without the planetary shadow, grain charges quickly converge to equilibrium values, 
and the stability map in Fig. [T3k looks very similar to one for a constant (positive) charge 
(Fig. [TJo). In Fig. [T3k . the superimposed bold-faced curves, corresponding to a +2.75 Volt 
constant potential, match the data very closely on far right. Since the electric potentials 
of the large grains rapidly converge to equilibrium, the Keple r-regime side of the radial 



i nstab ility also closely conforms to the analytical boundaries of 



Jontof-Hutter and Hamilton 



(I2OI2I ). The smaller grains, however, take significant amounts of time to reach charge 
equilibrium. The grains just to the left of the left-most solid curves in Fig. [T3k either 
escape (outside Rgyn) or fall into the planet (inside Rsyn), before they have enough time 
to reach their equilibrium charge. While these tiny grains experience modest electric 
charges, a different set of stability curves to the left of those in Fig. [T3k applies. Within 
1.2 Rp, for example, 0.01 fim sized grains collide with Jupiter within a few hours while the 
characteristic charging time is a day. Similarly, outside synchronous orbit, 0.01 /im-sized 
grains just outside the radially unstable zone for constant +2.75 Volt grains can still escape 
the planet. These grains, initially neutral, charge up slowly in the sunlight. Hence, even if 
their equilibrium charge would permit stable motion, the time spent in the radially unstable 
regime causes them to collide with Jupiter or escape before reaching charge equilibrium. 

With the planetary shadow turned on, the shadow resonance acts to increase 
eccentricities, destabilizing grains over a far broader range of sizes than than those that 
remain at their equilibrium potential. For the largest of these grains (a^ > Ifim), the 
shadow transit destabilizes grains over timescales commensurate with the precession due 
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Fig. 13. — Stability of grains with variable charge response to plasma for two grain charging 
models: a) constant photoelectron emission, and b) a more realistic photoelectric response 
which is interrupted during shadow transit. We adopt an aligned dipole for Jupiter's mag- 
netic field, and neglect its 3.12° obliquity. Grains, all launched at local noon, begin with zero 
electric potential and are integrated over 1 year, a timescale long enough to cover the orbital 
precession in the Kepler regime due to J2, which is included in our model of the gravity field. 
Consistent with the greyscale in Figs. |3] and HI the darkest region denotes radial instability- 
either collision with the planet or escape within Am = 5° of the equator plane- the moderate 
grey represents high-latitude collisions with Jupiter, and the lightest grey marks surviving 
grains with high latitude oscillations. The two filled squares correspond to the trajectories 
illustrated in Fig. [121 Note that the horizontal axis, now marking increasing grain radius as 
is no longer constant, spans twice the range as in previous figures. In both a) and b), the 
numerical boundaries are included for an assumed electric potential of +2.75 Volts. 
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to the gravitational J2 term (iHoranyi and BurnallQQll ). Thus we include the J2 term in 
our model for Fig. [121 and adjust launch speeds to ensure launch from a circular orbit. 
The smallest grains in Fig. [T^ respond to changes in the charging environment over a 
longer timescale. Thus, grains that survive the initial charging process (the stable grains 
of Fig. [T3k) reach an electric potential that deviates little from its mean over the orbital 
period, and hence the shadow has little effect on grains smaller than 0.01/im in size. This 
region of Fig. [T3b essentially matches Fig. [T]d in the Lorentz limit, with either vertical 
instability or stable high-latitude oscillations between 1.29 Rp and 1.70 Rp. 

Grains between O.Ol/xm and 0.1/im launched outside 1.2 Rp, but within Rsyni experience 
charge variations that cause them to spend some fraction of each orbit in a radially unstable 
regime. Eventually they strike the planet, although the timing for this is unpredictable. 
As we saw in Fig. Wlh . such grains experience random walks in radial location but do not 
cross synchronous orbit. Roughly half the grains in this region of Fig. [T3b collided with the 
planet at high latitudes. 

Grains larger than 0.5/im launched outside Rgyn in Fig. [T5b experience excited radial 



motions and ver t ical motion close to the radial sta bility boundary of Fig. [T3k (L* 



Hamilton 



1993a 



Jontof-Hutter and 



from Jupiter ( iHamilton and Kriigei 



Hamilton 



2' 



2OI2I ). which extends the Thebe ring away 



20081 ). Inside Rsyn-, Fig. shows that the shadow 
resonance destabilizes grains more than 10 times bigger than the largest grains destabilized 
with the shadow switched off. The boundary between stable and unstable here is determined 
not by the time of the integration but by the precession timescale due to the higher-order J2 
gravity field component at Jupiter {k. 0.25 years in the main ring, and longer further out). 

We emphasize that resonant charge variation on dust grains due to the photoelectric 
current clearly has an important effect on grain dynamics. Epicyclic motion also provokes 
resonant charge variations due to both radial gradients in the plasma properties and the 
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varying dust-plasma speed with gyrophase. Of the effects studied so far: launch speeds, 
realistic magnetic fields and time-variable electric charges, the latter appears to be the 
most important. Non-zero launch impulses, by contrast, are a minor effect on grain-orbit 
stability. The relative importance of the different effects, however, will vary dramatically 
with plasma properties. 

For example, in a dense plasma like the lo plasma torus, the equilibrium charge is 
always negative (as in Figs. [7^ and|8^) even in full sunlight. Spatial or temporal gradients 
in plasma properties can also have a profound effect. Since all of the parameters for 
charging are very uncertain, we leave a thorough study of these effects fo r another paper. 



The constant charge maps here and in 



Jontof-Hutter and Hamilton! (120121 ) are still relevant 



though, serving as a measure of the minimum instability in a given system. It is also a 
good approximation for planets with significant obliquities, where dusty rings spend much 
of their host's long orbital period in direct and uninterrupted sunlight. 

Accordingly we turn to the other major magnetospheres of the Solar System and 
construct stability maps for dust grains with constant charge in complex multipolar 
magnetic fields. 



5. Other planets 



After our detailed investigation of Jupiter, we are now in a position to map and 
interpret stability results for each of the magnetized planets in the Solar System. We begin 
with the planet with the simplest magnetic field, Saturn. 
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5.1. Saturn 



S aturn's full magn etic field can be described by an aligned dipole with a slight vertical 



201 



offset (iCao et al 
(730 = 0.0274 Gauss ( iConnerney et al. 



). We model the full field with gio = 0.2154, ^20 = 0.0164 and 



19841 ). Figure [T^ shows the stability map for 



negatively-charged grains, with the numerically-determined stability boundaries for an 
aligned and centered dipole included for comparison, as in Fig. HI 

As at Jupiter, (see Fig. H]) Saturn's dipole offset increases the instability of grains to 
vertical perturbations. This eliminates the stable zone close to the planet that we see for 
the aligned dipole case, and moves the outer vertical stability boundary significantly further 
from the planet. The effect is stronger at Saturn than at Jupiter due to its relative large (720 
term and to the larger Rp/ Rgyn at Saturn, making the planet a bigger target. By contrast, 
at Jupiter (see Fig. H]), a large locally-stable region in the Lorentz limit close to the surface 
and at high survives the inclusion of (720- 

For positively-charged grains, in Fig. [T3)d, the offset dipolar field causes the vertical 
instability to join the radial instability, as in Fig. Hb. However, a tiny island of globally 



stable grains survives near (L* 



3, rrJRp = 1.3). The radia 



stability boundaries for 



2OI2I ). match Saturn's full 



an aligned dipole field for Saturn (jJontof-Hutter and Hamilton 
magnetic field remarkably well. 

At Saturn we also see a slightly wider range of charge-to-mass ratios excited by the 2:1 
resonance between epicyclic and vertical motions, when compared to Jupiter (Fig. H^). This 
is due to the range of launch distances extending further out in units of Rsyn in Fig. [T4h . 
The accuracy of the theoretical curve matching this resonance vindicates the use of an 
aligned dipole approximation for Saturn's magnetic field to calculate radial and vertical 
orbital frequencies. 
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Fig. 14. — Stability of a) negative and b) positive Kepler-launched grains, followed for 0.1 
years in Saturn's full magnetic field. The greyscale matches that of Figs. |3]and|H The darkest 
region denotes grains that are radially unstable to escape or strike the planet (positive charges 
only), the moderate grey scale marks trajectories that were vertically unstable to climb out 
of the ring plane and strike the planet at high latitude, the light grey scale marks grains that 
are vertically unstable in the equatorial plane but remain globally bound with mirror points 
further than |Am| = 5° from the equator, and the white area marks locally stable orbits. 
Superimposed on the data are the numerically-determ ined stability boundaries for an alig ned 



2OI2I), as 



and centered dipolar magnetic field model for Saturn fiJontof-Hutter and Hamiltonl l! 
well as the curve marking the analytical 2:1 resonance between epicyclic and vertical motion, 
which closely tracks band of grains that reach high latitudes but remain bound. The white 
circles mark the local stability threshold at Rgyn, where = 2 ± y^- 
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The transition from grains that are lost to the vertical instability to those that remain 
bound in the B ring is at 1.70-Rp or ^ 102,000 km in the Lorentz limit (Fig. UM . This is 
close to a large increase in optical depth in the B ring that begins around 1. 72 Rp, and losses 



to ero sion may play a role in ring evolution across this boundary. 



Northrop and Connerney 



( 119871 ) argued for a link between the inner edge of B ring and the vertical stability boundary. 
Their model for vertical motion predicted all highly- charged grains to be unstable within 
1.54Rp at Saturn, close to the sharp inner edge of the B ring. Their model, however, 
did not provide a detailed mechanism and also invoked poorly understood electrostatic 



effects ( iNorthrop and Hill 



1983 



Northrop and Connerney 



19871 ). Voyager 2 data revealed 



another transition in optical depth between 1.63 Rp and 1.65 Rp (98,000 - 99,000 km), a few 
thousand kilometers inside the vertical stability boundary at Saturn with its full magnetic 
field configuration. While the proximity of this transition to the vertical stability boundary 
is intriguing, a detailed model to explain this congruency remains elusive. 



5.2. Earth 



The Earth's magnetic field is dominated by a dipole tilted by a moderate 11.4° from 
the axis of rotation. For our full-field m o dels, we use magnetic field coefficients out to 



octupole order from 



Roberts and SowardI (119721 ). For the Earth, qiq < and the magnetic 



field is inverted compared to all of the giant planets. Thus for the Earth, L^. > for 
negatively-charged grains, and it is these negative grains that suffer the radial instability 



(IJontof-Hutter and Hamilton 



2012). 



Figure [15] compares the stability of grains in Earth's full magnetic field to an aligned 
dipolar model. With an aligned dipolar field, the vertical instability at Earth in the Lorentz 
limit is local, leading to a region of high-latitude globally-stable oscillations. For the 
positive grains {L^ < 0, Fig. [T5b) this locally vertically unstable region curves towards the 
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planet as |L^,| decreases, and only a small range of grains launched near the surface between 

^ —0.1 and —1.0 collides with the planet. This changes very little with the inclusion 
of Earth's higher-order magnetic field terms, as Fig. [TS] indicates. More dramatically, 
the higher-order terms (primarily gu) expand the vertical instability in the Lorentz limit 
to further distances, almost to Rsyn- In this region, grains do in fact collide with the 
planet. The feature is very similar to what we saw at Jupiter in Figs. Eland El except that 
the global instability region is far narrower in the case of the Earth. Just as at Jupiter, 
tilting the magnetic field does not significantly move the inner vertical stability boundary 
of Fig. [TSb- For the negative grains. Fig. [T^ shows that the vertical instability is also 
displaced towards Rgyn and curves slightly upwards to merge with the region of radial 
instability. The expansion of the vertical instability for both positive and negative charges 



nearly to synchronous orbit has a unique benefi t in assisting ; 



debris from this crowded regio n of Earth orbit (iHoranyi et al. 



1997 



Valk and Lemaitre 



20081). 



le re r noval of dusty space 



1988 



Juhasz and Horanyi 



The radial instability in Fig. [T5k . however, looks very different at Earth than 
at Jupiter (compare Figs. [6l and [T5b) . Nevertheless, outside synchronous orbit, the 
higher-order magnetic field terms h ave little effect and the radial stability boundaries of 
Jontof-Hutter and Hamiltonl (120121 ) match the data remarkably well. Inside synchronous 
orbit, the more complicated magnetic field slightly extends the area of HRLOs that abut the 
disjoint regions of equatorially-confined radial instability. Thus Earth's full magnetic field 
barely alters the radial instabilities expected for an aligned dipole: The few unstable grains 
that exceed the radial stability boundary on the right side of Fig. [T5h may be associated 
with Lorentz resonances. In particular, the small cluster of points near (L^. = 0.01, 
^h/Rp = 2.0) closely corresponds to the inner 6:1 Lorentz resonance. All else being equal, 
the Lorentz resonances are more important close to the planet where the magnetic field 
irregularities are strongest. 
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Fig. 15. — Stability of Kepler launched negative (a) and positive (b) grains in Earth's 
magnetic field complete to octupole order, integrated over 1 year, with launch at 00 = 0. 
The solid curves mark the numeri cal stability boundaries for the Ear th with its anti-aligned 



Qio component alone (Fig. 12 from iJontof-Hutter and Hamilton 



20121 ). Grains in the darkest 



region, (negative charges only) crashed into the planet or escaped at a latitude less than twice 
Earth's 11.4° tilt angle. The moderate grey region marks grains that struck the planet at 
latitudes higher than twice the tilt angle, the light grey marks grains with stable vertical 
oscillations with mirror points exceeding twice the tilt angle in latitude, and white regions 
are the remaining stable orbits. The white circles mark the points = 2±\/3, at = Rsyn, 
where grains are on the threshold of radial instability. 



-48 - 



As an exercise, we compared the results displayed in Fig. [15] with a simpler tilted 
dipole model, including just the Qiq and gu terms (figure not shown). The main difference 
that arises is that the positive and negative grains that are excited to high latitudes, near 
= 10, tl/ Rp = 1.3) are not excited in the tilted dipole model. A more subtle difference 
is the extra set of collisions to the right of the radial instability boundary on the right-hand 
side of Fig. [T5b marking grains that were lost because of Lorentz resonances. These grains 
survive in the simple tilted dipole model. Deviations at greater distances are not expected 
due to the steep radial dependence of the quadrupole and octupole terms, and indeed, they 
are not seen. All in all, a tilted magnetic dipole is a robust model for the motion of charged 
dust grains at the Earth. 



5.3. Uranus 

Uranus' complex magnetic field destabilizes grains for a much wider range of 
charge-to-mass ratios than Jupiter, Saturn or Earth. Figure [16] shows the stability of grains 
launched a t Uran us, with magnetic field coefficients out to octupole order taken from 



Ness et al. 



(1l99ll ). This figure highlights the significant dependence of launch azimuth on 
grain lifetimes. We determined grain orbit stability for 12 equally-spaced azimuthal launch 
positions, and followed trajectories for 1 year. 

Beyond Rsyni the full Uranian magnetic field causes a large class of grains to escape 
rapidly, for both negatively- and positively-charged dust. In general the stability maps for 
negative and positive grains are very similar, especially inside synchronous orbit. Within 
Rsyni all trajectories in the Lorentz limit appear unstable for both negative and positive 
charges, which significantly constrains the low-energy plasma environment in the uranian 
ring system. Furthermore, both Figs. [T6k and [T6b show far more dependence on launch 
phase at high charge-to-mass ratio, on the left side of the stability maps, than on the right. 
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consistent with our results for Jupiter (Fig. [6]). 

Uranus' magnetic tilt as well as its quadrupole and octupole magnetic field coefficients 
are much more important when compared to its qiq term than at Jupiter or Saturn, causing 
strong Lorentz resonances, and hence a dramatic expansion of escaping negative grains 
over that seen at Jupiter in Fig. [TUb . Furthermore, at Uranus, grains as close as Rgyn 
can escape, unlike at Jupiter. In Fig. [T6k . a spike at (L* = — 0.03, ri/ Rp = 2) appears 
to be associated with the 2:1 inner Lorentz resonance that approaches the Kepler limit 
at 2.04 Rp. Interestingly, this resonance appears stronger for negative grains than for 
positive ones. At Jupiter, two spikes in Fig. [TUb distinguish the inner 2:1 vertical and radial 
resonances. Similarly, at Uranus, Fig. [T6b hints at an even stronger pairing of destabilizing 
2:1 resonances, one curving downward towards the planet as increases, and one arcing 
slightly upward. 



5.4. Neptune 



We model Neptune's magnetic field configuration with data from IConnerney et al. 



( 1l99ll ). As for Uranus, above, the stability map includes the effect of launch longitude on 
grain-orbit stability. And as with Fig. [161 Fig- [13 indicates the number of 12 equally-spaced 
launch azimuths that survive a 1-year integration. 

Figure [T7k maps the stability of negatively- charged dust at Neptune, and includes a 
large region of escaping negative grains, though this range of is slightly smaller at Neptune 
than at Uranus (Fig. [T6k). However, the escape region for negative grains is still much 
more significant at Neptune than at Jupiter (see Fig. [TOb). and it too reaches Rsyn- As 
at Uranus, grain orbit stability on the Lorentz-dominated side of Figs. [T7k and [T7b is 
strongly dependent on the launch phase. Inside synchronous orbit, grains in much of the 
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Fig. 16. — Stability of Kepler-launched grains in Uranus' full field over one year. There are 
three shades of grey, plus white to highlight the effect of azimuthal launch position on grain 
orbit stability. The darkest grey marks grains that were unstable for all of 12 equally-spaced 
launch longitudes. The intermediate grey denotes unstable trajectories for 6 to 11 launch 
positions, the lighter grey 1 to 5 launches, and the white regions were stable for all launch 
positions. The embedded curves mark the Uranian equivalent of those in Fig. [H the stability 
boundaries for an aligned dipole magnetic field model. The white circles mark the local 
radial stability at Rgyn, where = 2 ± -\/3. 
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highly-charged Lorentz hmit are unstable, except for a small region around ri/Rp = 2.2 
for both positive and negative grains, where stability varies significantly with launch 
phase. This contrasts with Uranus, where all grains inside Rgyn in the Lorentz regime were 
unstable. Again, however, the instability at Neptune vastly exceeds that of Jupiter. 

In the Kepler-dominated regime of Fig. [T7k . two spikes most likely associated with the 
inner 2:1 Lorentz resonance feature prominently. Just as we saw at Uranus, in Fig. [T7b. this 
inner 2:1 resonance appears to be weaker. For negative grains, the white stable zone to the 
right in Fig. [T7k reaches to higher values (smaller grain radii) than we saw for Uranus 
(Fig. [TBk). For both planets, the dependence on azimuthal launch position for stability is 
only important for \L^\ » 1. In the Kepler regime, grains move rapidly across magnetic 
field lines, and instabilities are effectively averaged over all launch phases. All evidence 
points to greater instability at Uranus than at Neptune. This is consistent with the Uranian 
dipole tilt of 59° exceeding Neptune's 47°. 



6. Discussion 

In this paper, we have studied three main effects on charged particle motion: i) 
non-zero launch velocities from orbiting parent bodies ii) complex magnetic fields and iii) 
time-variable electric charges. We presented data from over 250,000 numerical integrations 
and compared our results to analytical theories, extending the important concept of Lorentz 
Resonances to arbitrary charge-to-mass ratios, and show ed that effects of order Ll can 



explain results that the linear theory of iHamiltoru (119941 ) misses. 



Non-zero launch impulses (section 3.1), relative to the Kepler flow have little effect 
on charged-grain dynamics and stability. Radial stability boundaries are only noticeably 
affected by an azimuthal kick. Vertical instability, by contrast, is affected by vertical 
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Fig. 17. — Stability of Kepler-launched grains in Neptune's full magnetic field, for twelve 
equally-spaced azimuthal launch positions, and integrated over one year. The greyscale is 
as in Fig. [16] with the darkest grey unstable for all launch longitudes and grains in the white 
areas surviving for all launch positions. As in previous figures, we superpose numerically- 
determined stability results for an aligned dipole as solid curves, and mark the local radial 
stability threshold at Rgyn, where = 2 ± -\/3 with white circles. 
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impulses as might be expected, but also by azimuthal kicks which strengthen or weaken the 
magnetic mirror force. Finally, a radial kick barely affects dust grain motions at all. 

When considering the stability of grains in the wide variety of planetary magnetic fields 
in the Solar System, we have shown in section 4 that dust grains with constant charges 
provide the maximum stability possible at each planet. The simplest magnetic field that we 
have considered, that of Saturn, is well-described as an untilted dipole, moderately offset to 
the north. This offset noticeably expands the vertical instability but has little discernable 
effect on radial motions. 

Jupiter's magnetic field is substantially more complex than Saturn's with a moderate 
tilt, a southward offset, and sizeable higher-order field coefficients. The tilted dipole 
strongly affects vertical stability boundaries, and the loss of axisymmetry powers Lorentz 
Resonances. These resonances act to destabihze dust particles, allowing even negative 
grains outside synchronous orbit to escape from Jupiter. Thus the lo plasma torus, in 
which grains are expected to have negative charges, is not an impermeable barrier to 
escape. The high-speed dust streams detected by Ulysses and Gahleo near Jupiter likely 
originate in the lo torus. We extend Lorentz resonances from the Kepler regime by rewriting 
their frequencies in terms of the general radial, vertical, and azimuthal frequencies valid 
at all charge-to-mass ratios. This allows us to calculate the radial locations of Lorentz 
resonances as a function of L*. We note strong correlations between zones of instability 
and the predicted locations of Lorentz resonances. Finally, our numerical simulations show 
that some resonances with strengths proportional to Ll must be active. We show how to 
determine the frequencies and rough stengths of these high-order resonances. 

Our results for Jupiter are directly applicable to the Earth, which also has a magnetic 
field that is dominated by a moderately-tilted dipole. Two interesting differences, however, 
distinguish charged-particle motion at Earth from that at Jupiter. Firstly, due to its 
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inverted magnetic dipole, the radial instability at Earth affects negative not positive 
charges. Similarly, positively-charged dust at Earth behaves as negatively- charged dust at 
Jupiter. Secondly, the Earth is very small compared to the size of its synchronous orbital 
distance, minimizing the tendency of the vertical instability to force grains to collide with 
Earth. Furthermore, due to the rapid decay with distance of higher-order magnetic field 
terms, the effect of Earth's full magnetic field on dust grain trajectories differs little from 
that of a simple tilted dipole. 

Uranus and Neptune both have complex magnetic field configurations which render 
aligned or even simple tilted dipolar models insufficient. Both of these planets have 
substantial quadrupolar and octupolar components, which act to destabihze both negative 
and positive grains across the synchronous orbital distance, and over a far greater range of 
charge-to-mass ratios than at the other planets that we have studied. These distant planets 
highlight how increases in magnetic field complexity dramatically exacerbate dynamical 
instabilities. Future spacecraft missions will provide more detailed planetary magnetic field 
configurations than we have available today, but changes to the stability maps that we have 
provided here for const ant- charge dust grains are likely to be modest. 

Relaxing the assumption of constant charge-to-mass ratios also leads to a substantial 
increase in the range of dust-grain sizes that are destabilized. Results are highly model- 
dependent, and for simplicity we adopted a sparse plasma with constant spatial density 
and photoelectric charging. The time-variable charging currents on a dust grain due to 
passage through the planetary shadow significantly expand the range of grain sizes that 
are globally unstable, particularly inside synchronous orbit. Larger dust grains respond 
rapidly to changes in the charging environment and hence stability is determined by 
comp aring the destabilizing timescale of variable Lorentz fo rces with the orbital precession 



time (IHoranyi and Burns 



1991 



Hamilton and Kriiger 



20081 ). For our nominal model, this 
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increases the threshold radius for destabihzed grains by more than an order of magnitude. 
For smaller grains, by contrast, charging is slow with the result that different stability 
curves apply at different times, expanding the zone of instability by an order of magnitude 
over that expected for a constant charge. The sparse plasma that we have adopted here 
is appropriate for the dusty main ring and gossamer rings at Jupiter. We find that the 
removal of dust at Jupiter is dominated by the basic dipolar radial instability for positive 
grains, substantially extended to both larger and smaller particles by the effects of variable 
charging. We leave a full study of the dependence of these effects on plasma parameters for 
a future study. 
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